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We investigate and contrast the non-perturbative infra red structure of 
N = 1 and N = 2 supersymmetric non-compact U(l) gauge field theory in 
three space-time dimensions with N matter flavours. We study the Dyson- 
Schwinger equations in a general gauge using superfield formalism; this en- 
sures that supersymmetry is kept manifest, though leads to spurious infra red 
divergences which we have to avoid carefully. In the M = 1 case the superheld 
formalism allows us to choose a vertex which satisfies the U(l) Ward iden- 
tity exactly, and we find the expected critical behaviour in the wavefunction 
renormalization and strong evidence for the existence of a gauge indepen- 
dent dynamically generated mass, but with no evidence for a critical flavour 
number. We study the M = 2 model by dimensional reduction from four 
dimensional M = 1 electrodynamics, and we refine the old gauge dependence 
argument that there is no dynamical mass generation. We recognize that the 
refinement only holds after dimensional reduction. 



I. INTRODUCTION. 

Dynamical mass generation is an interesting strong-coupling phenomenon which requires 
non-perturbative methods, and has been studied extensively in A^-flavour non-compact 
U(l) gauge theory in three dimensions (QED3) In this non-supersymmetric model 

the dynamics drives the coupling to a non-trivial fixed point in the infra red 0,0]: this leads 
to dynamical mass generation (and hence chiral symmetry breaking) and critical behaviour 
in the normal (chirally symmetric) phase. 

The situation for supersymmetric versions of electrodynamics is less clear. In an early 
work on three dimensional M = 2 supersymmetric QED, it was argued using component 
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formalism that (at least in Landau gauge) a dynamical mass is generated |T[]. It was pointed 
out that the vanishing of the effective potential in supersymmetric field theories meant that 
the question of energetic favourability could not be answered simply; it was conjectured that 
the issue of the selection of the finite solution over the vanishing one might be resolved at 
the level of the effective action, and the Ward identities arising from the supersymmetry [|I| . 
However, it has been argued on the grounds of gauge dependence that a (non-perturbative) 
non-renormalization theorem in four dimensional A/" = 1 supersymmetric electrodynamics 
forbids the dynamical generation of mass [§. On dimensional reduction, four dimensional 
A/" = 1 supersymmetry reduces to three dimensional M = 2 extended supersymmetry, and 
so the non-renormalization theorem should also hold for the three dimensional extended 
model. Recently there has been speculation that the gauge dependence argument for the 
four dimensional model is not sound The analysis of the four dimensional model 

employs a simplified form of the U(l) Ward identity, and a heavily truncated vertex [§. 
This is because in four dimensions supersymmetry and U(l) gauge invariance require, in 
addition to the usual three point vertex, n-point vertices for all n > 3 (see section [VI 1| ). 
This problem does not arise in the three dimensional A/" = 1 model, where supersymmetry 
and U(l) gauge invariance require three- and four-point vertices only. Recently there has 
been evidence from numerical studies |TI| that the M = 2 model does not generate a mass 
dynamically, in line with the dimensional reduction of the non-renormalization theorem, 
but in contradistinction to the results of reference [^. 

An early work on AT = 1 supersymmetric QED in three dimensions demonstrated that 
troublesome logarithmic divergences present in the wavefunction renormalization of non- 
super symmetric QED3 are absent in the supersymmetric model, and it was argued that 
this allowed dynamical mass generation to take place for numbers of flavours less than a 
critical number |T2|, which was related to the critical flavour number of non-supersymmetric 
QED3 [0]. The analysis was done in component formalism, in Landau gauge and with a trivial 
three-point vertex (set to unity) which does not satisfy the U(l) Ward identity except in 
the limits of zero transferred momentum and trivial (unit) wavefunction renormalization. 

It is therefore of interest to look more closely at the infra red structure of both M = 1 
and M = 2 supersymmetric U(l) gauge fleld theory. The actual dynamically generated mass 
is the p — > limit of the mass function and is the pole of the matter two point correlation 
function; as such, it must be a gauge invariant object. Only if the p — limit leaves a 
gauge independent constant can it be said with certainty that a mass has been dynamically 
generated. It is therefore crucial to give a careful treatment of gauge dependence, and to 
look closely at the full vertex. Since supersymmetry cannot be dynamically broken in a 
supersymmetric U(l) gauge field theory |T^, we do not have to consider the possibility of 
different dynamical masses for each component of the multiplet. With these issues in mind, 
we choose to work in superfield formalism, which keeps supersymmetry manifest, and we 
will, as far as possible, work in a general gauge. In components computations (see e.g. 



reference |]Tj]) it is necessary to use the U(l) Ward identity and the Ward identities arising 
from supersymmetry to constrain the vertices: the advantage of the superfield formalism is 
that we only need to examine the U(l) Ward identity, for the formalism ensures that the 
supersymmetry Ward identities are satisfied. In particular, by using the superfield formalism 
we will be able to choose a three-point vertex which satisfies the U(l) Ward identity exactly 
in the A/" = 1 case. It is known [^,|T^ that superfield formalism for gauge theories results in 
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spurious infra red divergences arising from the propagation of a gauge artifact: by careful 
choice of approximations and, where necessary, gauge, we can avoid most of the infra red 
divergences and they do not spoil our results. 

We use the same approach to study the non-perturbative physics of both the J\f = 1 and 
M = 2 models: we construct a non-perturbative ansatz for the matter superfield propagator, 
and study the behaviour of the functions appearing therein through the Dyson-Schwinger 
equations. In principle there are Dyson-Schwinger equations for the matter superfields, the 
gauge superfield, and for the full vertex. The gauge superfield Dyson-Schwinger equation 
will be resummed to leading order in 1/N, where N is the number of matter flavours; this 
is exactly what is done in normal QED3. The vertex will be constrained by the U(l) Ward 
identity. When the results of the resummation and a choice for the vertex are put into the 
Dyson-Schwinger equation for the matter superfield, coupled integral equations for the self 
energy and wavefunction renormalization result, which have to be solved. As with normal 
QED3, approximations have to be sought in order to look for solutions of the equations. 

For the Af = 1 model we will present two approaches for studying the Dyson-Schwinger 
equations. First, we consider some of the simplest approximations, which are sufficient 
when applied to normal QED3 to demonstrate the existence of dynamical mass generation 
(at least in Landau gauge) 0]. The advantages of making these computationally convenient 
approximations are twofold: they allow us to work in a general gauge throughout, and admit 
conversion of the integral equation to an equivalent differential equation. The differential 
equation is simple to solve, and has a gauge invariant, constant solution in the p ^ 
limit. The second method incorporates the full vertex, consistent with the U(l) Ward 
identity. Ideally we would like to probe the gauge dependence in exactly the same way 
again; however, the richer structure of the integral equations arising from this approach only 
admit more direct methods. In both approaches we find the expected critical behaviour in 
the wavefunction renormalization [ITgHTHI 



Z{p)r^(£j\ 7~C?(l/iV), (1.1) 

and exhibit the possibility of a finite dynamically generated mass. In contrast to reference 



1^ we find no evidence for a critical flavour number, above which there is no dynamical 
mass generation. 

For the M = 2 model we again make some simple approximations which allow us to probe 
the gauge dependence. Again we find the expected critical behaviour in the wavefunction 
renormalization (|1.1|) , but this time we find that there is no gauge independent dynamically 
generated mass. This is based on a refinement of the gauge dependence argument of ref- 
erence P] ; we note that the refinement relies on the existence of the scale arising from the 
compactified dimension, so that the non-renormalization theorem for the four dimensional 
model may still be evaded P,p!0 



In section || we construct the action functional for a A/" = 1 supersymmetric theory of 
matter flavours interacting with a U(l) gauge field in superfield formalism. We also give the 
dressed propagators for the model which result from this action. We construct the U(l) Ward 
identity in section ^ using the BRS method, and in section ^ we construct the Dyson- 
Schwinger equations for the model, which we use to analyse the non-perturbative properties 
of the theory. We present and briefly discuss our simple computation in section |V|, in which 
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we are able to probe the gauge dependence of the mass function by making convenient 
approximations in the integral equations. Our more complete computation follows in section 
[VI| , where we choose a vertex which satisfies the U(l) Ward identity. Here we cannot use 
the elegant methods of section 0, but we are able to demonstrate that the solution obtained 
there persists in this more complete computation. To conclude this part of the paper, we also 
discuss the question of whether the finite mass solution is preferred to the vanishing mass 
solution, in the context of the suggestion that this can be answered by appealing to the 
effective action and the supersymmetry Ward identities. We find that no extra information 
can be obtained in this way, in contrast to the models of references |TU|-pT[] which, crucially, 
have extra constraints. 



We turn to the M = 2 extended model in section [VI1| , where we build an action functional 
and dressed propagators for four dimensional = 1 supersymmetric A^-flavour electrody- 
namics. We construct the Dyson-Schwinger equations for the model in section |VIII| and then 



compactify one dimension to obtain the coupled integral equations for the three dimensional 
M = 2 model we wish to analyse. We compute the wavefunction renormalization and refine 
the argument of reference |p for the non-existence of dynamical mass generation. 

We append discussion and concluding remarks in section and briefly contrast with 
the situation in non-supersymmetric QED3. In the appendices we collect some useful results 
and basic features of A/" = 1 supersymmetry in three and four dimensions, which are essential 
to the analysis of the models we consider. 

II. THE M = 1 ACTION. 

The M = 1 rigid superspace in three dimensions can be parameterized by the usual 
three space-time coordinates and two Grassmann-odd coordinates arranged in the Majorana 
spinor where a G {1,2} ||22|, which we denote collectively by the symbol z. We collect 



some basic features of this N = 1 superspace and superfields in appendix ^ Consider 
complex scalar "matter" superfields $, $* (mappings from superspace to C) which transform 
under local U(l) transformations in the familiar way: 

$(2) — ^ ^'{z') = e+^'^^(^) <l>{z), 

<l>*(z) — > $*'(2')= e~'^^^'\ (2.1) 

Here e is the (dimensionful) gauge coupling and G{z) must be a real scalar superfield (a 
mapping from superspace to R) to preserve the superfield nature of in order to build an 
action functional invariant under local U(l) transformations, we must construct covariant 
derivatives which transform in the same way as the matter superfields themselves: 

V„ — >V^ = e^"^(^)V„e-*^^("). (2.2) 

The covariant derivatives can be written in terms of a real spinor (superfield) connexion 
Ta as follows: 

Va = Da- ^e^«; (2.3) 
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here Da is the normal spinorial derivative, covariant with respect to supersymmetry trans- 



formations [22|. The connexion transforms in the usual way under infinitesimal local U(l) 



t r ansf or mat ions : 

Ta{z) T'^iz') = Ta{z) + DaG{z). (2.4) 

We wish to consider a model with Af = 1 supersymmetry, local U(l) gauge invariance 
and matter flavours. The required action then comprises three parts: the gauge invariant 
classical field strength term for the connexion r^, a (Lorentz) gauge fixing term and a locally 
U(l) invariant kinetic term for the matter superfields $ and $*: 



S = Sf^' + Sf + 5^; (2.5) 

^ciass _ f p (--d'^d^d'^dA r 



S^ = j d'xd'e (^-0 [V"<l>]* [V„<l>] . (2.6) 

We have included in the matter part an implicit sum over flavours, which do not interact 
with each other but interact with the same gauge field. The parameter ^ is the familiar gauge 
parameter; supersymmetric Feynman and Landau gauges are given by ^ = 1,0 respectively. 
From this action it is easy to derive the renormalized propagators for the matter fields and 
connexion: 

A„^(p; 12) = -^1^-1^ [{lH)Po.pD' - 5\12). (2.7) 

The gauge field propagator includes the effects of massless matter loops (to leading order in 
1/A^) in the vacuum polarization factor [|I| 

(i + «/H)-^ 

Both the spinor and scalar components of the matter superfield $ contribute to this cor- 
rection, and so the vacuum polarization is exactly twice that of non-supersymmetric U(l) 
gauge theory |1|, |I2| , and the effective ultra violet scale a is given by 



a = (2.8) 

The second rank antisymmetric symbol Cap appearing in the gauge field propagator acts 
as a metric for spinor indices [^, (see appendix ^) and we have written three-vectors in 
a convenient spinor notation: Pafi is a symmetric second rank spinor. In all the above an 
abbreviated superspace notation has been used where, for instance, 
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(2.9) 



The interaction piece of the action is given by: 

5iNT = j d?xSe (^-ic'"^ier„$D^<l>* + ^C'"^zer„$*L'/3$-^e2C'°%r^$*$y (2.10) 

The Feynman rules derived from this action are given at the end of appendix ^ Finally, 
note that we have not used any Wess-Zumino type gauge fixing for the connexion superfield; 
this is not a gauge invariant truncation (nor one which respects supersymmetry), and we 
wish, as far as is possible, to investigate the full gauge dependence of the non-perturbative 
correlation functions. 



III. BRS IN VARIANCE AND WARD IDENTITIES. 

Following the BRS approach, to compute the Ward identities for this model we promote 
the gauge fixed action ( p.5| ) to an action invariant under an enhanced gauge symmetry. The 
gauge fixing term in (|2.5| ) transforms under infinitesimal gauge transformations as follows: 

Jc^r, {D^D^Df") = -^G iUD^V^) . (3.1) 

To make a BRS invariant action, we add a ghost term to the action, whose transformation 
properties cancel those of the gauge fixing term: 

^ghost^ j ^3^^2q ^_unw), (3.2) 

where u, w are real Grassmann-odd scalar ghost superfields. The full action is then invariant 
under the extended gauge symmetry 

= ie (ijw) $, 
S^^* = -ie$* {i]w) , 

6r,Ta = Da {r]w) , 

Sr,w = 0; (3.3) 

here ?/ is a real Grassmann-odd number. Note that 5'^ = 0. The generating functional can 
now be written 

Z = j E)^ iD$* E>u E>w e'^+'^^-'^% (3.4) 

where 

Ss = j d^x d^9 {r^ + J$* + K'^T^ + au + rw) , 

A3 = y d^xd'^e {J*5r,^ + J5r,^* + K'^Sr^T^ + (j5r,U + T^^w) . (3.5) 
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The set {J, J* , K°',a,T} are superfield sources; the last three are Grassmann-odd valued 
and is a spinor. Using the fact that rj'^ = 0, the Ward identity immediately follows 



= / iD$ iD$* E>r^ E>u E>w e^^+^^= x 

j (fx £e (J* 5^$ + J(5^$* + K'^d^V^ + adr^u + t5^w)^ . (3.6) 

To interpret this in terms of correlation functions for the model, we construct the quantum 
effective action (the index c indicates a classical value): 

r r^, M„ w,] =-i\nZ [J, J\ ir", a, r] 

- / d^x (fe ( J$* + J*$c + K'^Tl + (TMe + ™c) ; (3.7) 



the sources are given by functional derivatives of the effective action 



" ' 

6T _ 5T _ 

r[$„$:,r^,M„w;J = jd'xd'y ($:(x)A-i(x-i/)$,(y) + r,(x)(A-i(x-y))"^r^(l/)) 
+ Jd'xd'yd'z (^^:ix)TUy)^Gix,y,z)D^<!>,iz) 

+ <|.,(x) r^(i/) |g(x, y, z)D^^:{z)) + ■■■ (3.9) 
Now we can carry out the functional integral in the identity ( ^.6] ) to obtain 

= ld'x (^e^$>. - ^e§-w,^. + D^^w, + ^^^^r^) . (3.10) 

Taking functional derivatives with respect to Wc, $c and we obtain the Ward identity in 
configuration space: 

= D^{y)G{x, y, z) + i/\-^{x - z)5{x - y) - iA"^(x - z)5{z - y), (3.11) 
which in momentum space reads 

D\p - q)Gip,p- q, q) = iA-\p) - tA'^q). (3.12) 
The advantage of using a superspace formulation is now apparent, for using the identity 

{D\p)y = -p', (3.13) 
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the Ward identity (|3.12| ) can be inverted. In fact, as we shall see, the structure of the 
Dyson-Schwinger equations is such that it can always be organized so that D^{p — q) acts 
to the left of the full vertex, and so the Ward identity can be used exactly as it appears in 



( p.l2| ). Of course, as with the non-supersymmetric U(l) theory, the Ward identity does not 



constrain the transverse part of the vertex, for which 

D\p-q)G{p,p-q,q) = Q. (3.14) 

Since, as discussed above, the integral equations can be arranged so that the full vertex only 
appears with an accompanying D^, this will not be a problem for our analysis. Moreover, 
in not having to invert the Ward identity, we will not introduce any extra kinematical 
singularities as p — *• g; this is known to cause problems in the non-supersymmetric model 



In principle the four-point vertex could also be constrained by the U(l) Ward identity; 
however, the form depends in a complicated way on the three-point vertex and it is sufficient 
for us to choose a trivial function for our four-point vertex. As we will show, the graph 
which includes the four-point vertex contributes only to the wavefunction renormalization, 
and the contribution vanishes in supersymmetric Feynman gauge. 



IV. THE DYSON SCHWINGER EQUATIONS. 

The non-perturbative properties of the model are determined through the Dyson- 
Schwinger equations. These equations lead to coupled integral equations for the self energy 
and wavefunction renormalization, which can be solved (in principle) to yield the dynami- 
cally generated mass: 

M = lim M(p) = lim ^^4- (4-1) 

The truncated Dyson-Schwinger equation we will use is shown graphically in figure |l|. 
The graph on the left hand side is a convenient shorthand for the difference between the full 
inverse propagator and the bare inverse propagator. 




p-q 




p p q 

FIG. 1. Schematic form of the Dyson-Schwinger equation for the non-perturbative matter 
two point correlation function. Solid lines represent matter superfield propagators, and wavy lines 
gauge superfield propagators; blobs indicate full non-perturbative quantities. 
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The graph on the left hand side of figure |l| is given by 



-i{Z{p) -I) j Sd 9)D^{p)^\p, 9) - zS(p) j So ^(-p, d). (4.2) 



The first "seagull" graph on the right is simply 



(2^ 



d^9id^92 <i>{-p,9i 



__G'4C^-52(12)A„^(P; 12) 



$*(P,^2), 



(4.3) 



and the last graph is given hj {P = p — q) 



d^q 
(2^ 



d^9i d^92 <^{-p,9i 



-Gs{p,P,q)D^{q)Df'{q)A{q-12) 



A^p{P;12)<i>*{p,92 



(4.4) 



Our use of a symmetric second rank spinor notation for three vectors immediately implies 
the following identities 



= 0, 

6'^p-q. 



(4.5) 



Note that the second of these identities indicates that the C symbol from the four-point 
vertex in the last graph projects out only the (1—0 component of the gauge field propagator, 
and indeed the entire graph would vanish in the supersymmetric Feynman gauge, ^ = 1. 

Upon substitution of the propagators (|2.7|), the superspace integrations can be performed 
using the identities above along with the relations ( |Al| ) and ( [A4|) in appendix 0. We 
will present two computations in what follows: first, a computation in which the simplest 
approximations are made: G^ip.P.q) is set to Z{q), and the bifurcation method is used, 
in which S is set to vanish in the denominators of all the kernels. These approximations 
are drastic but computationally extremely convenient, and are sufficient to demonstrate the 
existence of a dynamically generated mass in non-supersymmetric U(l) gauge theory [Q. 
These approximations will give us enough flexibility to probe the gauge dependence of the 
solution. Second, we present a more complete computation, in which the full three-vertex 
consistent with the Ward identity ( |3.12| ) is used. 

The treatments in the two cases follow a similar strategy: inserting the propagators ( p.7|) 
into equations (WM and ([4.3]) along with a choice for the three-point vertex, the superspace 



parts of the integrals can be unpacked and reassembled using the identities ( f4.5|) , ( |A1|) and 
to appear as functions of {Z, S,p, g} multiplying the two superspace structures: 



d^9 ^-p,9)D\p) ^*{p,9); 
d^9 ^-p,9)^*{p,9). 



(4.6) 
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Comparison with equation ( [4.2| ) shows that the function multiplying the first of these struc- 
tures is to be identified with the wavefunction renormalization, and that multiplying the 
second with the self energy function. From this point the usual procedure followed with the 
non-supersymmetric model can be adopted to study the coupled integral equations which 
multiply the structures above. The angular integrations can be performed, leaving integral 
equations now over the variable |g| only, which can in principle be evaluated. In practice, the 
kernels have to be approximated with low-momentum expansions. This is justified, for the 
(supersymmetric and non-supersymmetric) model is super renormalizable, and dynamical 
mass generation, if it occurs at all, occurs in the deep infra red. Two methods of solution can 
be employed for the resulting approximate integral equations: conversion to an equivalent 
differential equation, or direct integration with trial functions. The first of these will be used 
with our simple approximation, and the second for the more complete computation. 



V. SIMPLE COMPUTATION: GAUGE DEPENDENCE. 

In this section we present a computation of the wavefunction renormalization and self 
energy functions based on the simplest approximations. In defence of the approximations, 
they are sufficient to demonstrate the existence of a dynamically generated mass in non- 
supersymmetric U(l) gauge theory [Q; the three-point vertex we choose is compatible with 
the Ward identity only in the limit of vanishing transferred momentum (p), though it is 
consistent with a non-trivial wavefunction renormalization. The results of the computation 
for our model are clear: the expected critical behaviour |[I6|-[T8| is found for the wavefunction 
renormalization, and there is a simple form for a gauge independent dynamically generated 
mass. This is shown in a neat way, by converting the integral equation for the mass function 
to an equivalent differential equation, from which it is obvious that a ^-independent constant 
solution exists. It is impossible to use this neat method in the more complete computation 
to follow in section for the rich structure of the integral equation resulting from the full 
vertex prohibits the conversion to a differential equation. 

In the approximation G3 = Z{q) the superspace parts of the integrals can be unpacked 
as described above, yielding the following integral equations for the wavefunction renormal- 
ization and self energy function (K(g) = Z^{q) + S^(g)): 

^(^) - ^ / ' W / l^TT^2<^+«' ('"^ - • <=-^) 

where P"^ = + — 2pqx. The seagull graph ( |4.3| ) gives only an irrelevant p independent 
contribution to Z which would vanish in the supersymmetric Feynman gauge, and does not 
contribute at all to S: 

The angular x integrations in ( ^.1|) can be performed easily to obtain 
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Z{p) = 1 + 



1 a 



dq 2q^^^^^^^j^-^Ii{p, q; a), 
dq 2g2 



K(g) 



g; a) + (l-0X2(p, g; a)] 



where the kernels Xj are given by 



Xi = 

2pq 



p2 _ q2 
a 



+ ^ In 



\p - q\ p + q_ 

p + q + a 



p2 _ q2 



In 



p + q 
\p-q\ 



\p- q\+ a_ 



X2 = — In 

pq 



p + q + a 



]p- q\+ a 

In the hmit of p, g <^ a the logarithms can be expanded 
p + q + a 



In 



\p- q\+ a_ 



%-g)(^-2^ + ---)+(p-g), 



(5.3) 



(5.4) 



and upon using the bifurcation method in the denominators the gap equations take on the 
following approximate form: 



7( ^ 1 , 1 (1+0 r, 2 1(1 



+0 



S(p) 



(1-0 



c/g 



lS(g) 



pZ(g) A^7r2 7p g^(g) 



dq 



V g g^ 
lS(g) 



(5.5) 



A. Wavefunction Renormalization. 



The first of equations (|5.5|) can be integrated directly, to obtain the wavefunction renor- 
malization in the normal phase (where there is no mass generation) and as an approximation 
to the wavefunction renormalization when there is mass generation: 

^M-o +©' + /(!). (8-6) 

where a is a constant and the function / vanishes at least as quickly as p/a in the limit 
p — i> 0. In the above we have used the usual renormalization group argument to resum 
logarithmic terms as 

l + -,ln(^).( 

The exponent 7 is given by 



p\ 
a J 



7 -(1+0^. (5.7) 
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This critical behaviour is quahtatively the same as that found in non-supersymmetric elec- 
trodynamics [17,18], so this gives us confidence in our approach, and demonstrates that the 



rough qualitative behaviour is exhibited even within the drastic approximation scheme used 
here. The non-trivial wavefunction renormalization we have found here is to be contrasted 



with the early results of reference ||T2[ in which the wavefunction renormalization was taken 
to be unity. 



B. Self Energy. 

Turning now to the self energy function, we wish to construct an equivalent differential 
equation for the mass function 



this is achieved by differentiating the second integral equation (|5.5|) twice with respect to p 
and once with respect to ^. The reason we take an extra derivative with respect to ^ is that 
this will allow us to probe the gauge dependence of the solution directly, without having to 
fully solve the equation, for we are only interested in the gauge dependence of the solution. 
First we rescale variables a.s p ^ ap and M{p) i— > am(p) and then the resulting differential 
equation is: 

the superior point denotes a derivative with respect to C, and primes denote derivatives with 
respect to p. 

Demanding that m = for all values of p (in the infra red) yields the solution 

M(p) = Moe-P'/"', (5.9) 

which is obviously constant and ^-independent in the p ^ limit. We should be cautious, 
however, for the approximation of neglecting the mass function in the denominators might 
not be consistent with this limit. When, in the next section, we use the full vertex and 
restore the mass in the denominators, we will be able to show that this solution appears 
to persist in the limit p — >■ 0, confirming our assertion here. It is interesting to note that 
there are no solutions to (|5l8| ) which are gauge independent in the limit p ^ but gauge 
dependent elsewhere (this can be shown by developing a gauge dependent power series for 
m and then constraining the gauge dependence via the indicial equation): this contrasts the 
case of non-supersymmetric QED3, in which such solutions do appear to exist, and seem to 
lead to a critical fiavour number. We will return to a full discussion of this issue in section 
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C. Concluding Remarks. 



The approximations we have used in this section are rather drastic, but the results 
obtained are neat and transparent: the wavefunction renormahzation exhibits the expected 
critical behaviour [p!^-[l8 



7~^, (5.10) 

and we find a gauge independent solution for the dynamically generated mass 

M(p) =Moe-P'/"'. (5.11) 

These results are in line with those found in non-supersymmetric QED3, though we have no 
evidence here for a critical fiavour number above which no mass generation occurs, which is 
the case in the non-supersymmetric model. This contrasts the results of reference [12|, in 
which a critical fiavour number was found. We discuss this issue further in section ^ after 
we have presented our full computation. 

Had we not used the bifurcation method in this section, the resulting differential equation 
for the mass function would have been non-linear. What is curious is that the bifurcation 
method works so well: as we will show, the solutions derived above are entirely consistent 
with the full vertex and the restoration of the mass in the denominators of the kernel. 
The inaccuracies arising from regions of integration where the bifurcation method is least 
reliable are suppressed in the limits considered for the determination of the dynamical mass 
(i.e. vanishing p). 

VI. FULL COMPUTATION. 

In this section we turn to our more complete computation, where we use a three-point 
vertex consistent with the U(l) Ward identity ( |3.12D . The trick we use is to rewrite — 



—D^{P) in the gauge field propagator, and then integrate by parts until the D^{P) operators 
appear acting on the three-point vertex G3, so that the Ward identity (|3.12|) can be used 
as it appears. When this is done, the superspace parts of the integral equations can be 



unpacked as described in section IV and reconstructed after some algebra to obtain the 



coupled integral equations (in which we have already performed the angular integrations): 

^^^^ = ^ + iA^ / dq ^[^ilH)Zip)Ziq)J,ip,q;a) 

+ 2(1-0 {q'Z{p)Z{q) - q^Z\q) + S(p)S(g) - T?{q)) ^p, q; a)] , 

+ 2(1-0 (g'S(p)Z(g) -p2s(g)Z(p)) J2{p,q;a)] ; (6.1) 
the angular integrals Jli{p, q; a) are given by 
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2pq 
1 

pq 



a 



\p-q\ 



p^ — q^ — ^ 
+ - In 



a [|p — g| p + q 



p + q_ 
p + q + a 
_\p- q\ + a_ 

^ In 



\p-q\ 



p + q 
\p-q\ 



+ 



In 



p + q + a 
_\p- q\+ a_ 



(6.2) 



Again, the seagull graph gives an irrelevant p independent contribution to Z, equation ( ^.21) , 
and also vanishes in the supersymmetric Feynman gauge. Since much of our final analysis will 
have to be restricted to this gauge (or at least to approximations equivalent to considering 
this gauge), even if we were to make a p-dependent ansatz for G4, the contribution would 
vanish on account of the C"'^ tensor structure in the vertex projecting out only the (1 — 
component of the gauge field propagator. 

To derive the integral equations we have used the following Taylor expansion for operators 
of the form D'^{p + r) when acting on $*(p, 6): 



D\p) + -r-^ 



d 



J r=p 



1 

+ 2 



2 dr^pdr'^f^ 



(6.3) 



J r=p 



Note that this is an exact expansion, which terminates at second order. The curious extra 
factors of one half which appear in the expansion arise on account of the second of equations 
( [4.5| ) which puts a factor of two into the dot product of two three-vectors. When acting on 
$*(p, 6') we drop terms which involve derivatives with respect to r, for these correspond to 
higher derivative terms, which should be irrelevant for our study of the infra red physics. 

To study the integral equations it is again necessary to approximate the logarithms with 
small momentum expansions, and we obtain approximate integral equations for S and Z, 
in which we have used the abbreviated notations: 



A 



1 1 



,2 / 



Q{p, q; Z, S) = q'Z{q) [Z{p) - Z{q)] + E(g) [E(p) - S(g)] 

R{Z, T.) = i:{p)Z{q)-T.{q)Z{p) 
S{p, q- Z, S) = q^nv) Z{q) - p'^q), Z{p). 



(6.4) 



The wavef unction renormalization is given by 



«r Jo HQ) 



+ 4(1+0^ [ dq 
Jp 



Z{p)Z{q) f2q^ p' 
K(g) \ a a 



+2(1-0-4 r<;,«<^''^^'^>' ' 

Jo 



2g2 



pa. 

+2{l-i)A dq 

J V 



K(g) 
Q(p,g;Z,S) [ q^ 



K(g) 

and for the self energy we obtain 



p"^ — q^ ap"^ 



p q^ — p"^ 



(6.5) 
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The kernels are infra red singular; this arises from the propagation of a supersymmetric gauge 
artifact namely the lowest component of the gauge superfield; this singular behaviour 
in the infra red is a general problem for superspace formulations of gauge field theories, but 
is absent in supersymmetric Feynman gauge ^j. Note that the singular terms as p ^ g are 
regulated on account of the vanishing of the functions Q and S in this same limit. 

It is also evident from equations (|6.5| ) and (|6.6| ) that it is impossible to derive from them 



an equivalent differential equation (not even within the bifurcation method) on account of 
the presence of factors of Z{p) and in the kernels. Furthermore, it is evident that 
even if the factors of Z{p) and could be dealt with, the same problem arises from the 
{p^ — q^)^^ terms: the integrals cannot be entirely removed by differentiation. However, 
this last problem is removed in the supersymmetric Feynman gauge, in which the problem 
of the gauge artifact discussed above is eliminated: for this reason much of our analysis 
will be performed in Feynman gauge, or using approximations which are computationally 
equivalent to this gauge. 



A. Wavefunction Renormalization. 

In the normal phase, where S = 0, the integral equation for Z can be integrated directly. 
We begin with the approximation Z{j)) = Z{q), which is not quite as drastic as Z = 1, but 
has the same effect, namely to reduce the kernels to known functions of p and q. Performing 
the integrations yields the following form for Z (compare with equation ( |5.(j| )): 

^ 1 . ^ - m , (i±|) ,(1±|) 



aUIL)\f(l); (67) 



a/ \a 



where we have used the usual renormalization group argument to resum the logarithm, and 
as in equation ( ^^ , a is a constant of order 0{1/N) and the function / vanishes at least 
as fast as p/a in the limit p/a — >■ 0. The critical exponent 7 is found to have the usual 
behaviour. Note that the evaluation of the exponent 7 in this more reliable approach 
differs by a factor of two from the result obtained in section |V]; this is of little concern for 
the exponent is not a gauge invariant object. The dependence is the crucial property, 
and this is the same in both computations. 

We can demonstrate that this solution is stable, by now relaxing the assumption that 
Z{p) = Z{q), and feeding the solution 
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back into the integral equation. The parts of the integral equation which vanish in Feynman 
gauge now give divergent contributions, in the form of hypergeometric functions evaluated 
at the limit of their radius of convergence, and divergent logarithms from the remaining 
singularities at p = g. These singularities are again a consequence of the propagation 
of a supersymmetric gauge artifact; they can be successfully avoided by considering the 
supersymmetric Feynman gauge, ^ = 1, in which the iterated solution reads: 

7 = ^. (6.10) 

The critical exponent is unaltered, and the coefficients of the other terms receive corrections 
of order 0{1/N) in the denominator. Hence, for large N the solution is stable to iteration, 
at least in supersymmetric Feynman gauge. 



B. Self Energy. 



We turn now to the self energy, and the integral equation ( |6.6| ). As we discussed previ- 
ously, it is impossible to derive an equivalent differential equation, on account of the presence 
of and (p^ — q^)~^ terms in the kernels. We again study the equation in the Feynman 
gauge, which removes the difficulties of the singular terms. Dividing through by Z{p), we 
obtain the following integral equation for the mass function M{p): 



2A^7rV V7o Z{q) q"^ + M^{q) Z{q) q^ + M^{q) 

+ 1 M{p)-M{q) 2p' r 1 M{p)-M{q) 

^aj^ Z{q) q^ + M\q) a Z{q) q^ + M\q) ' ^ ' ' 

To study this we cannot use the usual method of setting M{p) = M{q) = M (0) everywhere 
and looking for self consistent solutions, for in this limit the kernels reduce to zero. Instead, 
we recall the solution of the differential equation in section 0, equation ( ^.9|) , and adopt the 
following ansatz for M inside the kernels: 

M(g) =moe-P'/"'. (6.12) 

We adopt also the solution 

for the wavefunction renormalization factors which appear above. Unfortunately the inte- 
grations cannot be performed analytically if the exponential factor in (|6.12|) is present in 



the denominators. Hence in the denominators we set M to a constant, m, and it is with this 
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in mind that the scale m will be interpreted as an average mass over the range p G [0,a]; 
we will see that in the limits we wish to consider, dependence on the introduced scale drops 
out. We can then perform the integrations analytically with result: 

11— ^\m) \\mj \ 1 1 \mj 

+l^e-^'/-' - ^ (^) e-V«^ esc ( '-^n) + f ^ 

21-7 4 \aJ V 2 J 21-7 Va 

1 7 /m\2^_p2/„2 _ 7 / 3-7 \ ^ / 7-1 ^ 

21 + 7Va/ 4Va/ V2yV2'a2 

^l^/my-y /ZV^''^^/. 3-7 5-7 fl^ 




23 — •y \ a J \ \ m J \ 2 2 \ m 



+4r(l^,(-)' )h-0(^), (6.13) 



where 



7 



The limit p — > is now well defined, and if we also assume m -C a and take the limit 7 <^ 1 
we find that 



M(0) 1 7 /vr 1 



mo 1 — 7 2V2 1 — 7 



+ ... (6.14) 



where the dots represent the omission of terms smaller than the finite terms exhibited above. 
It is clear then, that to leading order in 1/A^, the ansatz (|6.12|) is stable to iteration, and 
we have found that the solution of the differential equation in section |V| is consistent with 
the full vertex and the restoration of the mass in the denominators of the kernels. We have 
therefore demonstrated the possibility of a finite dynamically generated mass for this theory, 
though we still find no evidence for a critical fiavour number: we return to a discussion of 
this issue in section fX. 



C. Concluding Remarks. 

In this section we have presented a computation of the wavefunction renormalization 
and self energy with a three-point vertex which is consistent with the U(l) Ward identity. 
While we have not been able to use the elegant methods of section 0, we have been able to 
perform a direct computation, and within the approximations made we have found that the 
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wavefunction renormalization exhibits the expected critical behaviour |T6|-p!8[] and we have 
exhibited a finite solution to the gap equation for the mass function: 



M(p) =moe-P'/"'. (6.15) 

These results are qualitatively the same as those derived using the simple approximations of 
section 0, justifying the use of those approximations in determining qualitatively the gauge 
dependence in this model. 

Since we have found that a mass can be dynamically generated in this model without 
breaking supersymmetry, it is natural to ask whether it is possible to determine whether the 
massive or massless (for vtlq could still vanish) solution is favoured. In non-supersymmetric 
theories this question can be addressed by appealing to the effective potential (which upon 
using a suitable variational principle leads to the same non-perturbative physics for QED3 
as do the Dyson-Schwinger equations |]2^), but in supersymmetric models this potential 
vanishes. Based on experience with supersymmetric non-linear sigma models [|19| it was 



conjectured |I| that the question could be settled by appeahng to the effective action, and 
to the Ward identities arising from the supersymmetry. The same approach was taken 
in reference |2l| for a (2 + 1 dimensional) supersymmetric model with a CP^ constraint. 



The simplest Ward identity arising from supersymmetry relates the two-point correlation 
functions of the physical degrees of freedom in the matter superfield (in the presence of a 
mass m for the multiplet): 

(7/'„(x)?/'^(?/))o = (<9a/3 ^mCo^a) (0(a;)0(y))o- (6.16) 



Crucially, in the models of references [|l9|-pi[] , where there are further similar Ward identities 
relating correlation functions of fields arising from the implementation of the constraints, 
it was found that these Ward identities were only satisfied when the mass for the matter 
multiplet was non-vanishing. In the present model, and in the M = 2 model of reference [|1| , 
there are no such constraints, and so the Ward identities place no further requirements on 
the generated mass; of course, the superfield formalism we have adopted means that the 
Ward identity ( |6.16| ) is satisfied automatically. This leaves open the question of whether the 
vacuum selects the massless or massive solution in supersymmetric U(l) gauge field theory. 



VII. EXTENDED SUPERSYMMETRIC ACTION. 

The M = 2 extended rigid superspace is constructed in a way slightly different to the 
M = 1 model: it can be obtained by dimensional reduction from = 1 supersymmetry in 
four dimensions. The compactification is realised by considering all the fields in the model 
to be independent of one spatial dimension and then integrating out this variable. The 
rigid superspace of four dimensional M = 1 supersymmetry is parameterized by the usual 
four space-time coordinates and four Grassmann-odd coordinates arranged in two spinors 
9°" and where a G {1,2} and a G {1,2} p5|, which we again denote collectively by z. 
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This furnishes the superspace with two spinorial derivatives Da and -Dq,, covariant under 
super symmetry transformations, which leads to a crucial difference compared to the three 
dimensional model, in that the matter superfields now come in two types: chiral ($) and 
antichiral ($^), which respectively obey the constraints 



= 0, Da^^ = 0. 



(7.1) 



We collect some basic results for four dimensional superspace in appendix |B[ Consider chiral 
and antichiral superfields which transform under local U(l) transformations in the following 
way: 



(7.2) 



The functions A and A must respectively be chiral and antichiral, to preserve the chiral and 
antichiral nature of the matter superfields. Note that e is a dimensionless gauge coupling: 
the dimensionful parameter which will again lead to a dynamically generated scale will be 
the size of the compactified dimension. The gauge field is now a real scalar superfield, and 
transforms under infinitesimal U(l) transformations in the following way 



V{z) = V^iz) V\z') 



V{z) + - (A 



A). 



(7.3) 



We now wish to consider a model with M = \ supersymmetry, local U(l) gauge invariance 
and matter flavours; again the action functional comprises three parts, a locally U(l) 
invariant kinetic term for the matter superfields, a gauge invariant classical field strength 
for the gauge superfield, a (Lorentz) gauge fixing term; for convenience we will also add a 
bare mass term for the matter superfields: 



tGF 



I ok I on 
' m ' n 



(7.4) 



class 



s\ 



GF 



on 



(Tzd^d V 



1 



d^xd^B -V 
d^xd^e 
d'^xd'^9 fiQ $+$ 



V, 



d^x d^e i2o<!>\<^ 



t 



(7.5) 



Again we have included an implicit sum over flavours in the matter parts. There are a 
number of important differences between this action and the one we considered in the first 
part of this paper. First, there are now twice as many matter superfields: an antichiral 
and chiral superfield for each charge under U(l): as we shall see this changes considerably 
the structure of the propagators; in particular, note that in the absence of a bare mass 
//o the non-interacting correlation function ($_|_$_)o vanishes, and the Dyson-Schwinger 
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equations for the self energy and wavefunction renormalization are decoupled to an extent. 
Furthermore, it is evident that the implementation of the interaction with the gauge field is 
completely different: there are now an infinite number of vertices with increasing numbers of 
gauge superfields. In the Wess-Zumino gauge (which projects out the physical components 
of the gauge superfield, but which thereby breaks supersymmetry explicitly) the sequence 
of interaction vertices terminates at the four point vertex, giving it the same content as the 
model considered in the first part of this paper. In order to make the model tractable, we 
will truncate the Dyson-Schwinger equations at the level of the four-point vertex, effectively 
setting higher-point vertices to zero This situation is certainly not ideal [0, but we 
find that the graph in the Dyson-Schwinger equation arising from the four-point vertex 
contributes trivially to the wavefunction renormalization, and we therefore believe that this 
truncation will not affect the infra red physics. The dressed propagators for the model are 
as follows: 

A*±*±fo- 12) = ^^^^^ ^^^^ d^im (7 6) 

we have again used an abbreviated superspace notation where, for example 
5\12) = 6\e, - 02) = 5\e^ - 02) 5\h - h) = 5'(12) 5\V2). 



The undressed gauge superfield propagator reads |25 



A^^(j9; 12) = -Le^^^^"-^^*^^^ |^ (i _ ^) _ (1 + ^) 5\l2)^ ; (7.7) 

later we will dress this propagator with a vacuum polarization factor as before. 

For simplicity we will restrict our attention to the simple vertex approximation considered 
in the approach of section 0; this is in fact consistent with the U(l) Ward identity for this 
model p|, at least in the limit of vanishing transferred momentum p. 



VIII. THE DYSON SCHWINGER EQUATIONS. 

We can rewrite the matter superfields in terms of unconstrained scalar superfields, which 
makes the evaluation of the superspace parts of the Dyson-Schwinger equations simpler: 

^±{p,e,e) = e-'^'m±{p,e), 

^l{p^d,d) =e^^'^mi{p,d). (8.1) 

The Dyson-Schwinger equations for the current model are shown schematically in figures 
and ^. The graphs on the left hand side are a convenient abbreviation for the difference 
between the inverse full propagator and inverse bare propagator. 
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p-q 




P Q 



FIG. 2. Schematic form of the Dyson-Schwinger equation for the self energy function. Sohd 
hnes represent matter superfields, and the wavy hne represents the gauge superfield; blobs indicate 
full non-perturbative quantities. 




p-q 

p p q 

FIG. 3. Schematic form of the Dyson-Schwinger equation for the wave-function renormalization. 
Solid lines represent matter superfields and wavy lines represent gauge superfields; blobs indicate 
full non~perturbative quantities. 



The full self energy function on the left hand side of figure |^ is written as follows: 

J d^Oi d^92 ^+{-p,e2)j:ip)^5\i2)^4p,e^) 
d^e *+(-p,^)s(p)*_(p,^). 

The graph on the right of figure ^ is given by the following {P = p — q): 

d^q 



d'ed'e^+{-p,e2) [{-^e^G,) AH*l(g;i2)A^^(P;12)J ,<|._(j9,^i 
d'9d'9Ae'G,^^{-p,92) w^^?^Wt^^'(12) 



(27r 



2P2 



Z2(g) g2 + 4 

M;^(l-0-(l+05'(12)U-(l',^i)- 



(8.2) 



Expanding the exponentials and performing half of the superspace integrations using the 
identities in appendix O leaves: 



d^q 
(2^ 



4e^ r. 



_2p2 



:i-0 j d'e ^4-p,e)^4p,e). 



•4) 
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Note that in this model the superspace integrations project out the (1 — component of 
this graph. This is the basis of the gauge dependence argument |^ for the non-generation 
of dynamical mass in this theory. 

Turning to figure |], the graph on the left for the full wavefunction renormalization is 
written as follows: 
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d%d% e'^^'^^+i~p,92) (Zip) - 1) e-''^''-'^^''+^''^'' e'^^'^ ^^lipji) 
d% d% ^+i-p, 62) (Zip) - 1) e''^^'^ ^lip, e,) 

d'e {z{p) - 1) 6) e, e). (8.5) 



The superspace integrations can be performed in the same way as for the self energy to 
write the graphs on the right of figure |^. The "seagull" graph is computed to be a trivial p 
independent contribution, as in the model of the first part of this paper (and which again 
vanishes in supersymmetric Feynman gauge): 

^ d'e Ue'G, $+(-p, e, 6) A^^(g) ^{{p, 6, 6) 



d^q 



26^14-^(1-0 / d'9<!>iip,9,9)<^+i-p,9,9). 



.6) 



The last graph on the right of figure ^ is as follows: 



d\ 
(2^ 



462^3 



iZ{q) 



2P4 



-0-(i+0 



X / d'9<^>^{-p,9,9)<!>l{p,9^ 



(8.7) 



Collecting the results above, the Dyson-Schwinger equations of figures ^ and ^ then yield 
the following coupled integral equations for the self energy and wavefunction renormalization: 



Ziq) 1 

s^(g; 

jp + qf 

''Z\q) q^ + H^q) {p - qY 



(8.8) 



(8.9) 
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To study the infra red physics of A/" = 2 supcrsymmctric U(l) gauge field theory in three 
dimensions, we will simply consider all the fields above to be independent of one coordinate, 
and then integrate it out: this will introduce a natural scale p, the size of the compactified 
dimension, and we will then be able to compare directly with the model of the first part of 
this paper. We compactify one dimension as follows: 



/ 



1 

P 



(8.10) 



We can now incorporate the effects of massless matter loops into the gauge field propagator, 
by dressing it with the appropriate vacuum polarization factor: 



A^^(p;12) 



A^^ip; 12) 



l+a/\p\ ' 

where we have introduced the ultra violet scale 



(8.11) 



a — 



P 



in analogy with the treatment of normal QED3 and the A/" = 1 model considered in the first 
part of this paper. The combination e^/ p relates the dimensionless coupling e of the original 
model to the dimensionful coupling in the dimensionally reduced theory. We can carry out 
the angular part of the momentum integration easily in three dimensions, with result 
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p + q 
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p + q, 
p + q + a 
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(8.12) 



To study these equations it is once again necessary to expand the logarithms arising from 
the angular integrations, in the limit p,q <^ a, which leaves 



dq T{q) 



p2 _ q2 



q^ — p^ 



Zip) = 1+ 
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a 



j dqGiq) 
/ dqV, 



p2 _ q2 



ap 



e{p-q) + 



ap p"^ — q^ 



where we have introduced the convenient abbreviated notations 



q2 _ p2 



9{q-p) 



a 



Q{q) = G, 



Z\q) q^ + T.\q) 

Z{q) 
Z\q)q^ + T.\qr 



(8.13) 



(8.14) 



A. Wavefunction Renormalization. 

In the normal phase (or equivalently, within the bifurcation method) where S = 0, the 
integral equation for the wavefunction renormalization can be integrated directly upon using 
the simple approximation = Z{q), G4 = 1 M: 



Z{p) = 1+ 
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0? — p^ 
2pe 



(8.15) 



where the limit e — > 0"*" is implied. Again, outside the supersymmetric Feynman gauge, 
there are many infra red divergences which result from the propagation of a gauge artifact 
in the gauge superfield. However, ignoring these divergences as spurious, we can see that 
the solution above is consistent with the expected critical behaviour, with an exponent of 
order 0{l/N): 



Z{p) 



1 



2e 



(8.16) 



While it is not obvious that this form for the wavefunction renormalization is completely 
reliable, what is clear is that there is no evidence for an overall factor of (1 — the 
wavefunction renormalization, a fact which will be crucial for the discussion of the self 
energy which follows. 



B. Self Energy. 



It was argued in reference p[ that the overall factor of (1—^) in the first of equations (^.13|) 
and its four dimensional ancestor (|8.8|) for the self energy means that the self energy vanishes 
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in Feynman gauge, and since the full gauge dependence was carried through, this indicates 
that the self energy vanishes in all gauges; this was conjectured to be the result of a four 
dimensional (non-perturbative) non-renormalization theorem [Q. This argument has been 
criticized |P,|lO| on the grounds that the self energy function is not a gauge invariant function, 
and nor are the functions appearing in the kernel; in particular, the gauge dependence of 
the functions appearing in the kernel could in principle conspire to cancel the overall factor 
of (1—0 leave a gauge independent integral equation, which could have finite solutions. 
To answer this question it is convenient to divide the integral equation for S through by 
Z{p) to obtain 



M{p) 



a 



;i-0 / dq 



M{q) 



Z{p) g2 + M2(g) 



2q^p 

p2 _ g2 



e{p-q) + 



2p^q 
g2 _ jp. 



■9{q-p) 



(8.17) 



In the limit p ^ the function M{p) is gauge invariant; since we have shown that there 
can be no factor of (1 — from the wavefunction renormalization factor in ( |8.17| ), it seems 
that the function M{p), even if the result of integrating the kernel were finite, would retain 
the overall factor of (1-0; ^^id therefore the only self consistent solution would be that the 
mass vanished. We have not probed the structure of the vertex coming from the U(l) 
Ward identity [§], but this will only relate to combinations of S and Z: it will not lead 
to the gauge dependence required to cancel the overall factor of (1-0- In order to make 
this argument, we depend on the approximate form of the wavefunction renormalization 
( ^.l(j| ), and hence on the scale a; therefore our argument does not extend back to the four 
dimensional model where the scale a is absent, and the coupling remains dimensionless. 

A related point of interest is that, contrasting the A/" = 1 case, there is no part of the 
integral equation for M{p) which is not troubled by the infra red divergences generated by 
the gauge artifact; these divergences are known to disappear in supersymmetric Feynman 
gauge. The structure of the divergent terms make it impossible to construct an equivalent 
differential equation from the integral equation above. 



C. Concluding Remarks. 

In this part of the paper we have attempted to probe carefully the gauge dependence 
of the M = 2 model. We have chosen a simplified three-vertex and we have truncated 
the infinite series of higher-point vertices at four. Within these approximations, we have 
been able to compute the wavefunction renormalization, and have shown that although the 
computation is plagued with (expected) infra red divergences, we find no evidence for the 
overall gauge dependence which would be necessary to cancel the factor of (1-0 equation 
( ^.171 ) and evade the gauge dependence argument of reference [|]. We have followed this 
argument through for the (compactified) three dimensional model, where the scale a is 
available; in particular, the renormalization group resummed form for the wavefunction 
renormalization requires the existence of the scale a. While the properties of the four 
dimensional model are communicated to the three dimensional extended model through the 
compactification, it is not clear that the results presented here are applicable also for the 
four dimensional model where the scale a is necessarily absent; whence the arguments of 
references @,|T^ on evading the non-renormalization theorem in the four dimensional model 
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could still hold. Our results here are in agreement with a numerical evaluation of the effective 
potential for the three dimensional model in component formalism [11 . 



IX. DISCUSSION AND CONCLUSION. 

A. The J\f = 1 Model. 

The superfield formalism we have adopted has the advantage of keeping supersymmetry 
manifest; this has allowed us to concentrate on two aspects of the infra red physics which have 
troubled components computations in both the supersymmetric and non-supersymmetric 
versions of this model, namely gauge invariance and the full vertex. In particular, the super- 
field formalism gave us extra power to deal with the U(l) Ward identity for the full vertex. 
The disadvantages of the superfield approach are apparent in the spurious singular terms 
which appear as a result of propagating extra degrees of freedom in the gauge superfield. 
This problem can be avoided in supersymmetric Feynman gauge, and it is this which means 
we have had to give separate treatments to probe the gauge dependence and the issue of the 
full vertex. 

The first treatment employed some of the simplest approximations: a near trivial vertex 
(which satisfies the U(l) Ward identity only in the limit of vanishing transferred momen- 
tum, but which allows for a non-trivial wavefunction renormalization) and neglecting the 
mass in the denominators of the kernels. This programme of approximations is drastic, 
but computationally extremely convenient: enough to allow us to probe the gauge depen- 
dence of the mass function. This vertex choice and the bifurcation method have been used 
with some success in the study of non-supersymmetric QED3, and is sufficient to demon- 
strate the existence of a dynamically generated mass [Q] (at least in Landau gauge). In the 
model studied here, the method returns the expected critical behaviour in the wavefunction 
renormalization, and exhibits a gauge independent solution for the dynamically generated 
mass. The critical behaviour we find contrasts the trivial wavefunction renormalization of 
reference [jl2[, and is qualitatively the same as that of non-supersymmetric QED3, where 
mass is also dynamically generated. Our second treatment employs the biggest advantage of 
the superfield formalism, that of making the U(l) Ward identity tractable. Unfortunately, 
the problem of infra red divergences in the superspace approach restrict us to working in 
supersymmetric Feynman gauge for this computation. However, we find again the expected 
critical behaviour in the wavefunction renormalization and again exhibit the possibility of a 
finite dynamically generated mass; in particular, we find that the solution obtained in the 
first simplified approach is consistent with the full vertex, and a non-vanishing mass func- 
tion in the denominators. It is possible to criticize the approaches taken in both treatments 
we have presented, but since they tackle different issues with the same result, we argue 
that taken together they provide strong evidence for non-trivial infra red physics in A/" = 1 
supersymmetric U(l) gauge theory. 

The question of whether or not the finite mass solution is actually selected has been 
addressed: contrary to the suggestion in reference |l[] and the situations in the constrained 
models of references [|T9|- pl|] , we find that no extra information can be obtained by appeal- 
ing to the supersymmetry Ward identities. Unfortunately this leaves us with no way to 
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determine whether the finite mass solution is actually selected by the vacuum of this theory. 

We have found no evidence in this model for a critical fiavour number, above which 
no mass generation occurs. This is in contrast to the non-supersymmetric case 0,^. As 
has been noted already, the differential equation (|5.8| ) has finite solutions when the mass 
function is gauge invariant for all p (in the range permitted by our approximations: the infra 
red), but no solutions which are gauge invariant only in the limit of vanishing p and gauge 
dependent elsewhere. This is to be contrasted with the non-supersymmetric case, where 
solutions of the latter type do appear to exist. The integral equation for the mass function 
in that case (in general gauge, and with the simplified vertex choice of G3 = Z{q)) is as 
follows: 



M{p) 



1 q'M{q) r ( ia\ M{q) 



+ dq [q + 



p 2py q^ + M'^{q) J„ ^ y 2 J q^ + M^{q) 



(9.1) 



The conventional method for the non-supersymmetric case is to consider Landau gauge 
(^ = 0) and construct an equivalent differential equation |Q: 

d ( dM\ _ 8 p^M _ 

dp V dp J Ti^Np^ + M^' ^ ' ' 

the solution of this differential equation takes the form of a hypergeometric function; crucially 
this leads to a constraint on such that above a certain critical flavour number, no mass 
generation occurs 0]. Proceeding as we have done for the supersymmetric model, one has to 
make three differentiations with respect to p to remove the integrations, and the linearized 
(neglecting the mass in the denominators) result is as follows (in dimensionless variables 
rescaled by a): 



p^ m!" + m" + 



g 

6p + ^(p + 



Demanding as before that m = for all p, we flnd a constant solution for m; substituting a 
gauge dependent power series indicates that there is also the possibility of gauge independent 
constant as p — > but with a gauge dependent function of p away from this limit. This 
latter type of solution does not exist in the supersymmetric model. We see then that the 
gauge dependence of the two models is rather subtly different, and the fact that we flnd 
no critical flavour number may be related to this. In support of this view it is instructive 
to attempt to construct a differential equation similar to equation ( |9.2|) for the AA = 1 



supersymmetric model, which might then show some evidence for a critical flavour number. 
To this end we return to the second of integral equations ( ^.31) for the self energy, and 
again perform differentiations with respect to p, but retaining the mass as a constant in the 
denominators Q]: 

d_ ( , dM{p) \ ^ 2(1-40 v'M{p) 

dp V dp J Ntt^ p^ + M^' ^ ' ' 

By analogy with the analysis of reference 0] the "critical flavour number" derived from the 
equation above would be 
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8(1-40 



which on account of its exphcit gauge dependence has no physical meaning. 



An interesting comparison can be made with the model of reference [Q, where dynam- 
ical mass generation in an A^-flavour effective theory of fermions in three dimensions is 
studied using the renormalization group and the e expansion. The effective theory has for 
its degrees of freedom explicitly gauge independent (composite operator) fields, and so the 
treatment is necessarily gauge independent; in this approach also there is no critical flavour 
number. We also found no critical flavour number when we included the effects of the full 
vertex (and the non-trivial wavefunction renormalization): this is reminiscent of results in 



non-supersymmetric QED3 ||2^ where it has been argued that the critical flavour number of 
reference was an artifact of including only a trivial wavefunction renormalization and sim- 
plified vertex. In both our simple computation and the full computation, we have included 
the effects of non-trivial wavefunction renormalization. 



B. The J\f = 2 Model. 

In contrast to the M = 1 model, in extended supersymmetric electrodynamics, we find 
strong evidence for the non-existence of a dynamically generated mass, based on a refine- 
ment of the gauge dependence argument of reference [|]. We have used the same trunca- 
tion in the Dyson-Schwinger equations (in keeping only three- and four-point vertices) as 
used in the four dimensional model. This truncation is consistent with restriction to the 
Wess-Zumino gauge, which respects neither supersymmetry nor gauge invariance. Within 
this approximation we have chosen a simplified vertex (but one which is consistent with 
a non-trivial wavefunction renormalization) and we have then computed the wavefunction 
renormalization. We found no evidence for the appearance of an overall factor of (1 — 
which, if present, could have evaded 1^,0 the gauge dependence argument of reference [§. 



Instead we found that the wavefunction seems to have the same sort of critical behaviour as 
in the non-supersymmetric and M = 1 cases, though this is by no means certain on account 
of the presence of infra red divergences (again from an artifact in the gauge superfield). Our 



result here is in line with reference in which a numerical study showed that a mass 
is not dynamically generated in the M = 2 model. While we have been able to show that 
there is strong evidence against the dynamical generation of mass in the three dimensional 
extended supersymmetric model, our results depend crucially on the existence of the scale a 
(arising from the size of the compactified dimension). In the four dimensional model, where 
this scale is absent, the arguments of could still hold, and the non-renormalization 

theorem of reference HI be evaded. 
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APPENDIX A: SUPERFIELDS IN THREE DIMENSIONAL Af = 1 SUPERSPACE. 

Spinor indices are raised and lowered by the antisymmetric metric C^^, numerically equal 
to (72, such that for any spinor A: 

A^ = ^A^A^ ^ ]f^^ApA^. (Al) 

The spinorial derivatives, covariant with respect to supersymmetry, have the following 
explicit representation: 

Do.{p)^d^ + ePp^p, (A2) 

(where = ) and they obey the algebra 

{D^{p),Dp{q)} = {p + q)^p. (A3) 

The following identities follow directly from the definition of D: 

D'^iq) D^{q) = q'^^ + C^^D^q), 
6\12) (52(12) = 0, 
^2(12) (5^(12) = 0, 
^2(12) D^6\12) = (5^(12), 

D\p) S\12) D\q) 6\12) = D\p + q) 5\l2). (A4) 

Scalar superfields have the following component content: scalar, spinor and auxiliary 
fields, 

(t>{x)^<^{x,e)\ 

F{x) ^ D^^x,0)\ (A5) 

where the symbol | indicates evaluation at 6* = 0. Equivalently, the component content can 
be written: 

$(x, 6) = (l){x) + e^'ilJaix) - e'^F{x). (A6) 
The components of the connexion superfield are given by 
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U{x) 



r„(x,i 



^Df'D^r^{x,e)\. 



(A7) 



The components Ua and U can be gauged away in the so-called Wess-Zumino gauge. The re- 
maining physical components are the normal gauge field Vajs and its supersymmetric partner 
Afj. For further details, the interested reader is referred to reference p2[ . 



The superspace Feynman rules from the action in section ^ are as follows: 




■G,{p,p-q,q)C'^^Dp{q) 




-IGs{p,p-q,q) C^^Dpiq) 



P Q 




FIG. 4. Interaction vertices in = 1 model. 



APPENDIX B: SUPERFIELDS IN FOUR DIMENSIONAL M = 1 SUPERSPACE. 

Spinor indices are raised and lowered by the real metrics Ea/s and e^^, both numerically 
equal to the matrix —ia2, such that for spinor s x and r]: 

a a3 - -3 

X ^Xp Xa = e^px'^ 

XV = X'^Va XV = XaV""- (Bl) 

The spinorial derivatives have the explicit representation 

5a = -4-^^.4, (B2) 
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where {da, 0^} — 5^ and {da, 0^} — 5^. The covariant derivatives obey the algebra 

{Da{p),Da{q)} = {a'')^.{p + q),. (B3) 

The following identities follow directly from their definition: 

52(12)5^(12) = 6\12)6\12) = 
D\q)6\l2) = -4e(^i-^2)^e-i^2^j2), 

D\q)S\l2) = -Ae^'^^^'-^'^S^Oi - 62). (B4) 

Chiral superfields have the following component content (where now the symbol | now 
indicates evaluation at ^ = ^ = 0): 

(l){x) = ^x,e,e)\ 
^2Mx)-Da^x,e,e)\ 

F{x) = -^D'^x,9,9)\, (B5) 

which is equivalent to 

^x,e,e) = e-'^^^'^{x,e) 

^{x, 9) = (t){x) + ^29il]{x) + 9'^F{x). (B6) 

The component content for an antichiral superfield is similar, with 9 replaced by 9, and D 
replacing D. The component content of the gauge superfield is as follows: 

C{x) = V{x,9,9)\ 

Xa{x) = -iDaV(x,9,9)\ xc,{x) = DaV{x,9,9)\ 

{M{x) + iN{x)) = -D^V{x, 9, 9) \ {M{x) - iN{x)) = --D^V{x, 9, 9) \ 
\a{x) = -'-D^DaV{x,9,9)\ \a{x) = '-D''DaV{x,9,9)\ 
{a^)^.Vk{x) = -\[Da,D^] V{x,9,9)\ 
d{x)^\D''D''DaV{x,9,9)\ 

o 

45f d(x) - 2i {a^^aX dimVn] {x) = D^D^DpV{x, 9,9)\. (B7) 

In Wess-Zumino gauge the components C, x, X, ^ and N can all be gauged away, leaving 
the component expansion 

V{x, 9,9) = - {9a^9) Vk{x) + i9W\{x) - ¥9\{x) + -9Wd{x), 
V\x,9,9) = -\9'9W{x), 

V^{x,9,9) = Q. (B8) 
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For further details the interested reader is referred to [^. The Feynman rules are very 
simple; outside the Wess-Zumino gauge there are an infinite number of possible vertices, with 
two external matter superfield lines (one chiral and one antichiral), and n gauge superfields, 
the vertex corresponding to a factor of i{2e)"' (or i(— 2e)" for the fields $_ and $1., which 
are negatively charged under U(l)). In Wess-Zumino gauge, the series terminates at n = 2, 
and we dress the corresponding vertices with the functions G3 and G4. 
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